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We study the nonequilibrium spin transport through a quantum dot con- 
taining two spin levels coupled to the magnetic electrodes. A formula for the 
spin-dependent current is obtained and is applied to discuss the linear con- 
ductance and magnetoresistance in the interacting regime, where the so-called 
Kondo effect arises. We show that the Kondo resonance and the correlation- 
induced spin splitting of the dot levels may be systematically controlled by 
internal magnetization in the electrodes. As a result, when the electrodes are 
in parallel magnetic configuration, the linear conductance is characterized by 
two spin-resolved peaks. Furthermore, the presence of the spin-flip process in 
the dot splits the Kondo resonance into three peaks. 
PACS number(s): 73.23.-b, 73.63.-b, 75.25.+Z 
Key words: Kondo effect, magnetoresistance, quantum dot 
Spin-polarized transport in magnetic nanostructures, in particular, single-electron tun- 
neling in ferromagnetic (F) double junctions, has become a very active area of research, 
mainly because of its possible applications in information storage and processing devices 
|]J. In these junctions, the transport properties depend on the relative orientation of the 
magnetic moments of external electrodes. When the magnetic moments of the electrodes 
are antiparallel, the tunnel resistance increases. This is called as the tunnel magnetore- 
sistance (TMR). When the central grain in double junctions is small enough to form a 
quantum dot (QD), the effects of the discrete quantized energy levels as well as Coulomb 



blockade become significant, which has been considered theoretically and experimentally in 
the sequential tunneling regime P,p|,^j5|j6|J^,P,p|,^0|,p | . 



In the low temperature regime, a more subtle effect of large charging energies is the 



creation of new states of many-body character at the Fermi level by the Kondo effect [JT3 
For a QD coupled to the normal (N) electrodes, the Kondo effect is well understood in and 
out of equilibrium []l4, 15il6|, 17|. It is a consequence of a special kind of high-order tunneling 



process in which the electron inside the QD tunnels out followed by an electron with opposite 
spin tunneling into the QD. The whole system forms a spin singlet state and the net magnetic 
moment in the QD is zero. For the F-QD-F system, a very important question is, what are 
the consequence and characteristics of the Kondo effect when the magnetic moments of the 
electrodes are taken into account? 

In this Letter we have studied spin-dependent transport in an interacting QD coupled to 
two magnetic electrodes as shown schematically in Fig. 1. Different from the conventional 
Kondo problem in a N-QD-N system, here the characteristics induced by the strong elec- 
tronic correlation is sensitive to the relative orientation of magnetization between the two 
electrodes, namely, the parallel and antiparallel configurations as shown in Fig. 1. 

The model Hamiltonian for the F-QD-F system under consideration can be written as 

H = £dd„d a + Ud^d r d^d l + R(d+d l + H.c.) 

T a kaada + H.C. (1) 

kadL.Rcr kadL,Ra 

Here, the single particle energy e d is double degenerate in the spin index a, and the interac- 
tion is included through the Coulomb repulsion U. The first two terms in H represent the 
correlated level of the QD, the third term is used to describe potential spin-orbit coupling 
which may cause the spin rotation of an electron while in the QD. The spin-flip mechanisms 
that are relevant to the GaAs-based QD have been studied recently Jl2]| . The fourth term 
describes the free magnetic electrodes and the last term is the spin-dependent hybridiza- 
tion of the QD to the magnetic electrodes. Model (1) has been employed to study tunnel 
magnetoresistance in the sequential-tunneling regime [[Ll[ . 



Since the spin quantization axes in the electrodes are fixed by the internal magnetization 
of the magnets, the electrons tunnel into a superposition of spin-up and spin-down states. 
This coherence has to be taken into account when calculating current pofl . Technically, we 
introduce a spin rotation transformation = (l/y / 2)(c| =F Cj), in terms of which the dot 
Hamiltonian in Eq. (1) is rewritten as J2a £ caC^c a + Uc^c^c^ci with e ca = Ed ± R for up 
and down spins, respectively. The current through the left electrode can be calculated from 
the time evolution of the occupation number Nl = Y^ka ^XLa a kLa for electrons in the left 
electrode using nonequilibrium Green functions. The result is 



^ = e(iV L ) = |E/^Tr[V fci G< L (e)], 



(2) 



where we have defined the tunneling amplitude matrix 
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a = L, R, and lesser Green functions [G^ L {t)] aa i = i{a^ La (0)c a '(t)). Next, we use Dyson's 
equation to calculate the non-equilibrium Green functions and express J by the Green 
functions G r of the dot as follows 



J=IJ ^Tr{T L {e)[iG<{e) + &(e)A(e)]}, 



(4) 



where [G^(t)] CTCT , = i(c^,(t)c a (0)) is the matrix expression for the lesser Green functions. 
A(e) = i[G^(e) — G^(e)] is the spectral function. f a (e) is the Fermi-distribution function 
in the a electrode, and 



T a (e) 



r?(e) + r?(e) T?(e) 



(5) 



is the line-width matrix with T"(e) = 27rJ]fcga |Hicra-| 2 <5(e — Cfcao-)- The spin-dependence of 
T°(e) originates from the bulk magnetization of the electrodes. 

In order to determine the retarded Green functions we choose the equation of motion 
(EOM) method. Although this method does not provide quantitative results in the Kondo 



regime, it gives the qualitative feature both in equilibrium and nonequilibrium cases |T6 . 
The method generates higher-order Green functions, which have to be truncated to close 
the equation. In the infinite-?/ limit we obtain 

Gf (c) = [el -e c - E*(e) - Sf (^(I - n c ) (6) 

where I is a 2 x 2 unit matrix, {i c ) ua ' = 5 ua ie ca , T, R = —i(T L + T R )/2 is the self-energy 



matrix due to the tunneling coupling between the electrodes and dot, Sf is due to the many- 

r \Vka]\ 2 f(^ka]) _|_ |Vfca|| 2 /(<Sfcc.x) 1 



body correlation, with the matrix elements [Sf (e)] CT(T = E fca [ j_ yfcaTl _ /(efc " T) + l^aJ /(£fcai) 



and [£f(e)U = E fca [ ' Vfc "£ /(eteT) - |Vta £f tei) ] ; (n c W = <W(c+e<f) must be calculated 
self-consistently. 



To solve the lesser Green functions Gf we generalize Ng's ansatz |19[ to the present case. 
The interacting lesser and greater self-energies are assumed to be of the form S <!> = Sq ^B, 
where B is a matrix to be determined by the condition Ti < — S > = S R — S A . This ansatz 
is exact in the non-interacting limit (U = 0) and guarantees automatically the current 
conservation law. As a result one obtains S < = Sq [S^ — S^] _1 [S i? — S ]. Using this ansatz, 
Gf can be obtained by Keldysh equation Gf = Gf S < G^ 1 . Substituting the expressions 
of the QD's Green functions in Eq. (4), defining S < = T R (T L +T R y\i: R - and after 
a straightforward calculation, one obtains a compact expression of the tunneling current 

J = \ I |^[r L Gf ^G^e) - f R (e)} (7) 

This expression generalizes the current formula in Ref. [16] to the spin-dependent Anderson 
model with additional spin-flip relaxation and allows one to describe the coherent spin 
transport through an interacting quantum dot coupled to magnetic electrodes. 

In the following calculations, for simplicity we neglect the energy dependence in the tun- 
neling matrix elements. The intrinsic linewidth of the dot energies has a form T°(e) = 
T"8(W — |e|) with the electrode band width 2W 3> max(fceT, eV, T a ). We consider 
two magnetic configurations, namely, parallel and antiparallel configurations. When the 
magnetic electrodes are in parallel configuration, we assume that the spin-majority elec- 
trons are up (a =|) and the spin-minority electrons are down (a =!). We further assume 



that in the antiparallel configuration the magnetization of the right electrode is reversed. 
Therefore the spin dependence of the coherent transport can be conveniently considered 
by introducing magnetic polarization factors pl and pn for the left and right barriers, re- 
spectively. = r (l ±Pl), Tf(i) — a r (l ±Pr) is for the parallel configuration, and 
■T'fd) = ar (l =F Pr) for the antiparallel configuration. T describes the coupling between 
the quantum dot and the left electrode without internal magnetization and a denotes tun- 
nel asymmetry between the left and right barriers. In this work we assume the symmetric 
barriers, i.e., a — 1, pl — Pr — P- 

The spin-resolved spectral densities are calculated via the relation Pt(i) ( e ) = 
-±Im{(Gf) TT + (Gf) u =F (Gf) n =F (Gf) iT }. Figure 2 shows the behavior of the spec- 
tral densities in parallel and antiparallel magnetic configurations, which will be used to 
discuss the conductance results below. The Kondo resonance for each spin is clearly mani- 
fested by a sharp peak at e = (the chemical potential is set to be zero) in spectral densities 
for both magnetic configurations. However, the peak shape is sensitive to the magnetic 
configurations of the electrodes. As observed from Fig. 2(a), in parallel configuration the 
excitation characteristics of the spectral densities are remarkably different from the normal 
case in two prominent ways: (i) the Kondo resonance for down spin is enhanced, while the 
up-spin resonance is suppressed; (ii) interestingly, the broad single-particle resonances shift 
away from the dot level e^, but in opposite directions for different spins. The excitations for 
the down spin shifts towards higher energy, while it shifts to lower energy for the up spin. 
This splitting is due to the spin-dependence of the interacting self-energy matrix Sf(e), 
whose real part is different for the up and down spins and sensitive to the values of the spin 
polarization factor p. From the expression of Sf (e), the renormalized spin levels Sda are 
given by the self-consistent equation 

~ \ ^ l^fcairl fi^kaa) 
£da — £d + 2^ ~ ~ ~ 

k a ^dcr ^kaa 

2irk B T 1 e d -ii a 

= £d + ^ 2^ {ln -IT + Re * [ 2 - 1 2^T ]} ' (8) 

where \& is the digamma function. The result of Eq. (8) is shown in the inset of Fig. 2 (a), 
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where the dressed dot levels e&j are plotted as a function of p. This spin splitting of the 
dot levels in the interacting regime due to the magnetic properties of the electrodes leads 
to essential changes in the transport properties (see below). In addition, as p increases, 
the spectral weight of pj.(e) goes up, while the spectral weight of Pt(e) is reduced. Thus 
a net magnetic moment is induced by the magnetic coupling. In the case of antiparallel 
configuration, the spectral densities of two spins are identical and the Kondo peak is not 
influenced by the presence of magnetic polarization of the electrodes [see Fig. 2(b)]. 

The level dressing and Kondo resonances in Fig. 2(a) are suppressed when taking into 
account the spin-flip process. Remarkably, it shows in Fig. 3(a) that a large spin-flip 
transition R splits the original Kondo peak into three well-defined peaks. As seen from the 
expression of Sf (e), besides the peak at e = 0, the additional two peaks appear at e = ±i?, 
respectively. It is different from the antiparallel case with symmetric barriers, in which, as 
in the normal case, only the two Kondo peaks evolves from the presence of the spin-flip 
process [see Fig. 3(b)]. 

The dramatic changes in the spectral densities when rotating magnetic moments of the 
electrodes from parallel to antiparallel configuration suggests substantial different transport 
properties for these two configurations. Figure 4 shows the linear response conductance G as 
a function of Ed, which can be tuned via the external gate voltage, for different temperatures. 
In the antiparallel configuration, as observed from Fig. 4(a), the temperature dependence 
of the conductance is similar to the normal case although with a lower amplitude: the 
Kondo resonance broadens the conductance peak and saturates the peak amplitude at low 
temperatures. In addition, the peak shape remains nearly symmetric over a broad range 
of temperatures. In the parallel configuration, however, the conductance peak becomes 
asymmetric with decreasing temperature as shown in Fig. 4(b). This asymmetry is even 
more pronounced in the spin-resolved conductance G a [see the inset of Fig. 4(b)]. The 
peak splitting in G-\ and G± is in the same manner as shown in Fig. 2 due to the magnetic 
dressing of spin levels. Thus their superposition results in a double-peak structure at low 
temperatures and for large p as shown in Fig. 4(b) (doted line). The main peak with larger 
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amplitude is dominated by the up-spin resonance, while the other peak with lower amplitude 
nearly comes from the down-spin resonance. This novel spin filtering effect is fully caused 
by the interplay between the dot correlation and the magnetic coupling. When spin-flip 
transition is included, as shown in Fig. 4(c), the splitting of the conductance is suppressed. 

To describe the dramatic change of G when the magnetic moments in the electrodes are 
rotated from antiparallel to parallel configuration, we define the linear magnetoresistance as 
MR = (G p — Gap) /G a p, where G p and G ap are the linear conductance in the parallel and the 
antiparallel configurations, respectively. Figure 5 shows the magnetoresistance as a function 
of Ed for different temperatures. One can see that when the dot level is far from resonance, 
the magnetoresistances approach to the same value for different temperatures, coinciding 
with the non-interacting case. At low temperatures, however, dramatic changes occur in the 
resonant tunneling regime. The magnetoresistance develops into a dip with a negative value. 
This means that the linear magnetoresistance significantly decreases and even changes its 
sign at low gate voltages. In the empty orbital regime where the dot level is higher than the 
chemical potential of the electrodes, the magnetoresistance may be enhanced to a value as 
large as 160%. 

In summary, using Anderson model, it is shown that the magnetic moment arrangement 
in the electrodes plays an essential role in spin-dependent transport of a F-QD-F system in 
the interacting regime. For parallel magnetic configuration, the Kondo resonance and QD 
energy levels can be controlled by the magnetic polarization in the electrodes. Consequently, 
the linear conductance appears as a spin-resolved double-peak structure and a net magnetic 
moment emerges in the QD. The spin-flip process in the QD results in a splitting of the 
Kondo peak into three peaks. We expect these results are useful in exploiting the role of 
electronic correlation in spintronics. 

This work is supported by CNSF under Grant No. 69625608 and by US-DOE. 

Note added. — After the work was completed, we noticed that Kondo problem in TMR 
has also been considered by Sergueev et al. JH]] . However, they did not discuss spin splitting 
of the conductance and spin-flip effects, which compose essential points in this paper. 
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Figure captions 

Fig. I. Schematic plot of the F-QD-F system considered in this work. 

Fig. 2. P|(e) (solid line) and pj.(e) (dotted line) as a function of e in the (a) parallel, and 
(b) antiparallel magnetic configurations for p = 0.5, /c#T = 0.02r , e rf = — 4r , and R = 0. 
Inset in (a) indicates spin splitting of the dot levels as a function of p. 

Fig. 3. p(e) as a function of e in the (a) parallel, and (b) antiparallel magnetic configu- 
rations with R = 0.2r . Other parameters are the same as in Fig. 2. 

Fig. 4. Linear conductance G as a function of Ed in the (a) antiparallel, and (b) parallel 
magnetic configurations in the absence of spin- flip process with p = 0.5. Inset indicates 
spin-resolved conductance G^ and G^. The effect of of intradot spin-flipping on the parallel 
conductance is shown in (c). 

Fig. 5. Linear magnetoresistance as a function of for several values of temperature. 
Other parameters are the same as in Fig. 4. 
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Fig. 1 




Fig. 3 




Fig. 4 



